Search for: All records

Abstract Population III (Pop III) stars, the first generation of stars formed from primordial gas, played a fundamental role in shaping the early Universe through their influence on cosmic reionization, early chemical enrichment, and the formation of the first galaxies. However, to date, they have eluded direct detection due to their short lifetimes and high redshifts. The launch of the James Webb Space Telescope (JWST) has revolutionized observational capabilities, providing the opportunity to detect Pop III stars via caustic lensing, where strong gravitational lensing magnifies individual stars to observable levels. This prospect makes it compelling to develop accurate models for their spectral characteristics to distinguish them from other stellar populations. Previous studies have focused on computing the spectral properties of nonrotating, zero-age main-sequence (ZAMS) Pop III stars. In this work, we expand upon these efforts by incorporating the effects of stellar rotation and post-ZAMS evolution into spectral calculations. We use the JWST bands and magnitude limits to identify the optimal observing conditions, both for isolated stars, as well as for small star clusters. We find that, while rotation does not appreciably change the observability at ZAMS, the subsequent evolution can significantly brighten the stars, making the most massive ones potentially visible with only moderate lensing. 

Abstract The disks of active galactic nuclei (AGNs) are expected to be populated by numerous stars, either formed in the outer regions of the disk via gravitational instability or captured from the nearby nuclear star cluster. Regardless of their formation mechanism, these stars experience altered evolutionary paths, mostly shaped by the accretion of dense disk material. In this study, through the comparison of different timescales, we chart the evolutionary outcomes of these AGN stars as a function of disk radius and across a range of supermassive black hole masses, spanning from 10⁶to 10⁹M_⊙, for two popular AGN disk models. We find that in the outer regions of the disk, stars evolve similarly to those in the interstellar medium, but in the inner and denser regions, accretion quickly turns low-mass stars into massive stars, and their fate depends on just how quickly they accrete. If accretion occurs at a faster rate than nuclear burning, they can reach a quasi-steady “immortal” state. If stars accrete faster than they can thermally adjust, runaway accretion occurs, potentially preventing a quasi-steady state and altering the disk structure. During the AGN lifetime, in the regions of the disk that produce massive stars, supernovae (SNe) and gamma-ray bursts (GRBs) may occur within the disk over a wide range of optical depths and ambient densities. Subsequently, in the final phase of the AGN, as the disk becomes depleted, formerly immortal stars will be unable to replenish their fuel, leading to additional SNe and GRBs. 

Abstract Relativistic jets from a Kerr black hole (BH) following the core collapse of a massive star (“collapsar”) is a leading model for gamma-ray bursts (GRBs). However, the two key ingredients for a Blandford–Znajek-powered jet—rapid rotation and a strong magnetic field—seem mutually exclusive. Strong fields in the progenitor star’s core transport angular momentum outward more quickly, slowing down the core before collapse. Through innovative multidisciplinary modeling, we first use MESA stellar evolution models followed to core collapse to explicitly show that the small length scale of the instabilities—likely responsible for angular momentum transport in the core (e.g., Tayler–Spruit)—results in a lownetmagnetic flux fed to the BH horizon, far too small to power GRB jets. Instead, we propose a novel scenario in which collapsar BHs acquire their magnetic “hair” from their progenitor proto–neutron star (PNS), which is likely highly magnetized from an internal dynamo. We evaluate the conditions for the BH accretion disk to pin the PNS magnetosphere to its horizon immediately after the collapse. Our results show that the PNS spin-down energy released before collapse matches the kinetic energy of Type Ic-BL supernovae, while the nascent BH’s spin and magnetic flux produce jets consistent with observed GRB characteristics. We map our MESA models to 3D general-relativistic magnetohydrodynamic simulations and confirm that accretion disks confine the strong magnetic flux initiated near a rotating BH, enabling the launch of successful GRB jets, whereas a slower-spinning BH or one without a disk fails to do so. 

Abstract Although stable neutron stars (NSs) can in principle exist down to massesM_ns≈ 0.1M_⊙, standard models of stellar core-collapse predict a robust lower limitM_ns≳ 1.2M_⊙, roughly commensurate with the Chandrasekhar massM_Chof the progenitor’s iron core (electron fractionY_e≈ 0.5). However, this limit may be circumvented in sufficiently dense neutron-rich environments (Y_e< 0.5) for which $M_{\mathrm{Ch}}\propto Y_e^2$is reduced to ≲1M_⊙. Such physical conditions could arise in the black hole accretion disks formed from the collapse of rapidly rotating stars (“collapsars”), as a result of gravitational instabilities and cooling-induced fragmentation, similar to models for planet formation in protostellar disks. We confirm that the conditions to form subsolar-mass NS (ssNS) may be marginally satisfied in the outer regions of massive neutrino-cooled collapsar disks. If the disk fragments into multiple ssNSs, their subsequent coalescence offers a channel for precipitating subsolar mass LIGO/Virgo gravitational-wave mergers that does not implicate primordial black holes. The model makes several additional predictions: (1) ∼Hz frequency Doppler modulation of the ssNS-merger gravitational-wave signals due to the binary’s orbital motion in the disk; (2) at least one additional gravitational-wave event (coincident within ≲hours), from the coalescence of the ssNS-merger remnant(s) with the central black hole; (3) an associated gamma-ray burst and supernova counterpart, the latter boosted in energy and enriched withr-process elements from the NS merger(s) embedded within the exploding stellar envelope (“kilonovae inside a supernova”). 

Abstract Eruptive mass loss in massive stars is known to occur, but the mechanism(s) are not yet well understood. One proposed physical explanation appeals to opacity-driven super-Eddington luminosities in stellar envelopes. Here, we present a 1D model for eruptive mass loss and implement this model in theMESAstellar evolution code. The model identifies regions in the star where the energy associated with a local super-Eddington luminosity exceeds the binding energy of the overlaying envelope. The material above such regions is ejected from the star. Stars with initial masses of 10−100M_⊙at solar and SMC metallicities are evolved through core helium burning, with and without this new eruptive mass-loss scheme. We find that eruptive mass loss of up to ∼10⁻²M_⊙yr⁻¹can be driven by this mechanism, and occurs in a vertical band on the H-R diagram between $3.5\lesssim \log \left(T_{\mathrm{eff}}/\mathrm{K}\right)\lesssim 4.0$. This predicted eruptive mass loss prevents stars of initial masses ≳20M_⊙from evolving to become red supergiants (RSGs), with the stars instead ending their lives as blue supergiants, and offers a possible explanation for the observed lack of RSGs in that mass regime. 

Abstract Nuclear star clusters (NSCs), made up of a dense concentration of stars and the compact objects they leave behind, are ubiquitous in the central regions of galaxies surrounding the central supermassive black hole (SMBH). Close interactions between stars and stellar-mass black holes (sBHs) lead to tidal disruption events (TDEs). We uncover an interesting new phenomenon: for a subset of these, the unbound debris (to the sBH) remains bound to the SMBH, accreting at a later time, thus giving rise to a second flare. We compute the rate of such events and find them ranging within 10⁻⁶–10⁻³yr⁻¹gal⁻¹for SMBH mass ≃10⁶–10⁹M_⊙. Time delays between the two flares spread over a wide range, from less than a year to hundreds of years. The temporal evolution of the light curves of the second flare can vary between the standardt^−5/3power law to much steeper decays, providing a natural explanation for observed light curves in tension with the classical TDE model. Our predictions have implications for learning about NSC properties and calibrating its sBH population. Some double flares may be electromagnetic counterparts to LISA extreme-mass-ratio inspiral sources. Another important implication is the possible existence of TDE-like events in very massive SMBHs, where TDEs are not expected. Such flares can affect spin measurements relying on TDEs in the upper SMBH range. 

ABSTRACT Stars embedded in active galactic nucleus (AGN) discs or captured by them may scatter onto the supermassive black hole (SMBH), leading to a tidal disruption event (TDE). Using the moving-mesh hydrodynamics simulations with arepo, we investigate the dependence of debris properties in in-plane TDEs in AGN discs on the disc density and the orientation of stellar orbits relative to the disc gas (pro- and retro-grade). Key findings are: (1) Debris experiences continuous perturbations from the disc gas, which can result in significant and continuous changes in debris energy and angular momentum compared to ‘naked’ TDEs. (2) Above a critical density of a disc around an SMBH with mass M• [ρcrit ∼ 10−8 g cm−3 (M•/106 M⊙)−2.5] for retrograde stars, both bound and unbound debris is fully mixed into the disc. The density threshold for no bound debris return, inhibiting the accretion component of TDEs, is $$\rho _{\rm crit,bound} \sim 10^{-9}{\rm g~cm^{-3}}(M_{\bullet }/10^{6}\, {\rm M}_{\odot })^{-2.5}$$. (3) Observationally, AGN-TDEs transition from resembling naked TDEs in the limit of ρdisc ≲ 10−2ρcrit,bound to fully muffled TDEs with associated inner disc state changes at ρdisc ≳ ρcrit,bound, with a superposition of AGN + TDE in between. Stellar or remnant passages themselves can significantly perturb the inner disc. This can lead to an immediate X-ray signature and optically detectable inner disc state changes, potentially contributing to the changing-look AGN phenomenon. (4) Debris mixing can enrich the average disc metallicity over time if the star’s metallicity exceeds that of the disc gas. We point out that signatures of AGN-TDEs may be found in large AGN surveys. 

Abstract Observations show an almost ubiquitous presence of extra mixing in low-mass upper giant branch stars. The most commonly invoked explanation for this is thermohaline mixing. One-dimensional stellar evolution models include various prescriptions for thermohaline mixing, but the use of observational data directly to discriminate between thermohaline prescriptions has thus far been limited. Here, we propose a new framework to facilitate direct comparison: using carbon-to-nitrogen measurements from the Sloan Digital Sky Survey-IV APOGEE survey as a probe of mixing and a fluid parameter known as the reduced density ratio from one-dimensional stellar evolution programs, we compare the observed amount of extra mixing on the upper giant branch to predicted trends from three-dimensional fluid dynamics simulations. Using this method, we are able to empirically constrain how mixing efficiency should vary with the reduced density ratio. We find the observed amount of extra mixing is strongly correlated with the reduced density ratio and that trends between reduced density ratio and fundamental stellar parameters are robust across choices for modeling prescription. We show that stars with available mixing data tend to have relatively low density ratios, which should inform the regimes selected for future simulation efforts. Finally, we show that there is increased mixing at low reduced density ratios, which is consistent with current hydrodynamical models of thermohaline mixing. The introduction of this framework sets a new standard for theoretical modeling efforts, as validation for not only the amount of extra mixing, but trends between the degree of extra mixing and fundamental stellar parameters is now possible. 

Abstract Planetary engulfment events have long been proposed as a lithium (Li) enrichment mechanism contributing to the population of Li-rich giants ( A (Li) ≥ 1.5 dex). Using MESA stellar models and A (Li) abundance measurements obtained by the GALAH survey, we calculate the strength and observability of the surface Li enrichment signature produced by the engulfment of a hot Jupiter (HJ). We consider solar-metallicity stars in the mass range of 1–2 M ⊙ and the Li supplied by a HJ of 1.0 M J . We explore engulfment events that occur near the main-sequence turn-off (MSTO) and out to orbital separations of R ⋆ ∼ 0.1 au = 22 R ⊙ . We map our results onto the Hertzsprung–Russell Diagram, revealing the statistical significance and survival time of Li enrichment. We identify the parameter space of masses and evolutionary phases where the engulfment of a HJ can lead to Li enrichment signatures at a 5 σ confidence level and with meteoritic abundance strengths. The most compelling strengths and survival times of engulfment-derived Li enrichment are found among host stars of 1.4 M ⊙ near the MSTO. Our calculations indicate that planetary engulfment is not a viable enrichment pathway for stars that have evolved beyond the subgiant branch. For these sources, observed Li enhancements are likely to be produced by other mechanisms, such as the Cameron–Fowler process or the accretion of material from an asymptotic giant branch companion. Our results do not account for second-order effects, such as extra mixing processes, which can further dilute Li enrichment signatures. 

Abstract Stars are likely embedded in the gas disks of active galactic nuclei (AGN). Theoretical models predict that in the inner regions of the disk, these stars accrete rapidly, with fresh gas replenishing hydrogen in their cores faster than it is burned into helium, effectively stalling their evolution at hydrogen burning. We produce order-of-magnitude estimates of the number of such stars in a fiducial AGN disk. We find numbers of order 10^2–4, confined to the innerr_cap∼ 3000r_s∼ 0.03 pc. These stars can profoundly alter the chemistry of AGN disks, enriching them in helium and depleting them in hydrogen, both by order-unity amounts. We further consider mergers between these stars and other disk objects, suggesting that star–star mergers result in rapid mass loss from the remnant to restore an equilibrium mass, while star–compact object mergers may result in exotic outcomes and even host binary black hole mergers within themselves. Finally, we examine how these stars react as the disk dissipates toward the end of its life, and find that they may return mass to the disk fast enough to extend its lifetime by a factor of several and/or may drive powerful outflows from the disk. Post-AGN, these stars rapidly lose mass and form a population of stellar mass black holes around 10M_⊙. Due to the complex and uncertain interactions between embedded stars and the disk, their plausible ubiquity, and their order-unity impact on disk structure and evolution, they must be included in realistic disk models.